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We propose trial wave functions for quasiparticle and exciton excitations of the Moore- Read Pfaffian 
fractional quantum Hall states, both for bosons and for fermions, and study these numerically 
Our construction of trial wave functions employs a picture of the bosonic Moore-Read state as a 
symmetrized double layer composite fermion state. We obtain the number of independent angular 
momentum multiplets of quasiparticle and exciton trial states for systems of up to 20 electrons. We 
find that the counting for quasielectrons at large angular momentum on the sphere matches that 
expected from the CFT which describes the Moore-Read state's boundary theory. In particular, 
the counting for quasielectrons is the same as for quasiholes, in accordance with the idea that the 
CFT describing both sides of the FQH plateau should be the same. We also show that our trial 
wave functions have good overlaps with exact wave functions obtained using various interactions, 
including second Landau level Coulomb interactions and the 3-body delta interaction for which the 
Pfaflfian states and their quasiholes are exact ground states. We discuss how these results relate 
to recent work by Sreejith et al. on a similar set of trial wave functions for excitations over the 
Pfaffian state as well as to earlier work by Hansson et al., which has produced trial wave functions 
for quasiparticles based on conformal field theory methods and by Bernevig and Haldane, which 
produced trial wave functions based on clustering properties and 'squeezing'. 



I. INTRODUCTION 

The fractional quantum Hall plateau observed^ at 
filling = I has recently been at the center of much 
excitement, because it is expected that the elementary 
charged excitations of this state may be non-Abelian 
anyons.^"'' Moreover, manipulation of these anyons po- 
tentially represents an avenue to topologically fault toler- 
ant and hence scalable quantum computation^^ ^. Var- 
ious experiments have already observed important sig- 
natures of these excitations, such as their fractionalized 
charge^^ and tantalizing interferometric properties'^' 
which may provide a smoking gun detection of non- 
Abelian anyonic statistics'^ 

Moore and Read's Pfaffian wave function^, and its 
particle-hole conjugate, the anti-Pfaflftan'^'^^ remain the 
leading candidates for the description of the electronic 
ground state at = |. However, since Hall plateaus 
are probed through their excitations, it is of great im- 
portance to understand not only the ground state, but 
also the spectrum of low lying excitations which are nat- 
urally associated with Moore and Read's proposal. Much 
eflFort has been devoted to understanding the quasihole 
excitations, that is, the low energy states which are ap- 
propriate in situations where there is more magnetic flux 
piercing the sample than at the center of the Hall plateau 
(or equivalently, the electron density is lowered compared 
to the center of the plateau). A natural set of candidate 
wave functions exists for these, namely the exact zero 
energy ground state wave functions of the model 3-body 
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Hamiltonian introduced in Refs. 4,21. This Hamiltonian 
has the property that the Pfaflaan state is its highest 
density (or lowest angular momentum) zero energy state. 
At lower densities (when quasiholes are present), it has 
multiple zero energy states and one may conjecture that 
these states are good trial wave functions for the low 
energy Coulomb spectrum. In fact, the trial wave func- 
tions for localized quasiholes for the MR-states, which 
are constructed using correlators in conformal field the- 
ory (CFT), are automatically zero energy eigenstates of 
the 3-body Hamiltonian^, and hence coherent superpo- 
sitions of its zero modes, because this Hamiltonian en- 
codes properties of the operator product expansion in 
the CFT. Numerical studies^'^^'^^ on small systems have 
indeed found that the zero modes of the 3-body interac- 
tion provide a reasonable description of the low energy 
spectrum. 

Much progress has also been made in gaining a analyti- 
cal understanding of the zero energy states of this model 
Hamiltonian and its generalizations with /c-body inter- 
actions, which play a similar role for the Read-Rezayi 
series of states^^. Notably, it is known exactly how many 
independent zero energy states exist at any number of 
particles N and flux quanta 7V(/> and even how many of 
these states exist with any given angular momentum^^'^^. 
These countings are an important fingerprint of statisti- 
cal properties of the excitations and of the conformal field 
theory (CFT) which describes the edge of a system with 
boundary in the thermodynamic limit. 

Not nearly as much is known about the quasielectron 
excitations, which occur at higher density or lower mag- 
netic field, and about the neutral excitations (excitons), 
which can be viewed as combinations of quasiholes and 
quasielectrons. Clearly one may still conjecture (as done 
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implicitly e.g. in Ref. 6) that the low energy states of 
the 3-body Hamiltonian are good trial wave functions, 
but these are now no longer exact zero modes, and no 
exact expression is known for them. The aim of this pa- 
per is to propose and study an alternative set of trial 
wave functions for the neutral and charged excitations 
of the fermionic and bosonic Moore- Read Pfaffian states, 
for which we can write down an explicit analytical form. 
The fermionic wave functions we propose are very simi- 
lar in construction to those recently studied by Sreejith, 
T6ke,W6js and Jain in Refs. 27,28, which appeared while 
the present work was being written. In fact, though there 
are subtle differences between our construction and theirs 
(see Section IV) , which can cause differences in the qual- 
ity of the approximation of real systems, we believe the 
two sets of wave functions should be able to describe the 
same fractional quantum Hall universality class. Never- 
theless, our results on the counting of multiplets of quasi- 
electrons are markedly different from those presented in 
Ref. 28. There it was stated that the counting for quasi- 
electrons is different from that for quasiholes and the pos- 
sibility was raised that the quasielectrons might have dif- 
ferent braiding properties from quasiholes. However, we 
find that quasielectrons satisfy the same universal count- 
ing properties as quasiholes, suggesting that they also 
enjoy the same braiding properties. 

Our construction of trial wave functions is based on 
the idea that we can view the bosonic Moore- Read state 
as a double layer^^ of bosonic Laughlin 1/2 states — the 
fermionic MR ground state can be obtained by multi- 
plying with a Jastrow factor involving all the particles. 
A construction of localized quasiholes from the double 
layer representation of the MR wave function was al- 
ready given in Ref. 29. Here, we generate trial wave 
functions for all excitations of the bosonic and fermionic 
MR-states by constructing, in each layer, all the possi- 
ble excitations (excitons, quasiholes, and quasielectrons) 
over the z/ = 1/2 Laughlin ground states. The latter are 
well-understood in terms of the composite fermion (CF) 
theory that is, by creating quasielectrons and/or quasi- 
hole excitations in the integer quantum Hall effect of the 
CF at effective filling fraction = 1. We calculate the 
numbers of independent angular momentum multiplets 
of such states which exist at numerically accessible N 
and A^0, providing evidence that our construction yields 
quasielectrons with the same type of non-Abelian statis- 
tics and edge CFT as the quasiholes. We also compare 
the trial wave functions to the low energy eigenfunctions 
of the 3-body Hamiltonian and the second Landau level 
Coulomb Hamiltonian, showing directly that they indeed 
are good candidates to describe the excitations of the MR 
Pfaffian state. 

Other candidates for quasielectron and/or exciton trial 
wave functions, not directly based on composite bosons 
or fermions, have been proposed previously by Hansson, 
Hermanns, Regnault and Viefers^^'^^ and by Bernevig 
and Haldane^^. The construction in Refs. 31,32 is based 
on CFT and looks superficially very different to ours, but 



nevertheless, we expect that the quasiparticle wave func- 
tions presented there are the maximally localized (coher- 
ent) states, which can be produced from the states we 
propose here. The construction in Ref. 33 determines 
the trial wave functions for quasielectrons by requiring 
that they vanish when certain patterns of clusters of elec- 
trons are formed and also that they be dominated by cer- 
tain root configurations. We review both constructions 
in some detail in the rest of the paper and comment on 
the similarities and differences to our own construction. 



We have focused our numerical tests on systems with 
an even number of electrons, in part because these sys- 
tems exhibit a unique incompressible ground state, mak- 
ing quasielectrons and excitons clearly defined. Recently 
a number of works has also appeared which study the 
band of low energy states which appears in systems at 

= I when the number of electrons is odd^^'^^'^^. In 
particular, the paper by Sreejith et al.^^ employs trial 
wave functions for the states in this band which are based 
on a double layer CF system. Further recent work which 
focuses on the properties of excitons at z/ = |, and par- 
ticularly on the roton minimum, includes Refs. 36,37. 



Outline of the paper. In Section II, we give an 
overview of the Moore-Read state and of the quasielec- 
tron constructions of Hansson et al. and Bernevig and 
Haldane. In Section HI, we give a quick review of the con- 
struction of excitations over composite fermion or com- 
posite boson ground states. In Section IV, we describe 
our own trial wave functions for quasielectrons and exci- 
tons and explain how they can be numerically evaluated 
and studied both by real space Monte Carlo methods and 
by Fock space methods (using eigenstates of angular mo- 
mentum), which allow us to work at machine precision. 
In particular, we give details on the calculation of over- 
laps and of the number of independent trial wave func- 
tions for given N^jy and number of excitons or quasi- 
electrons. In this section, we also give a detailed explana- 
tion of the relation between our trial wave functions and 
those of Ref. 28. In Section V, we present our numeri- 
cal results, which include state counting of quasielectrons 
and excitons and overlaps between our sets of wave func- 
tions for bosons and fermions and the low lying states in 
the exact spectra obtained for bosons and fermions with 
three-body hardcore interactions, two-body hardcore in- 
teractions and second Landau level Coulomb interactions 
(with a slight shift of the pseudopotential Vi to obtain a 
stable Pfaffian state) . We note that our results for quasi- 
electron counting are consistent with the idea that the 
conformal field theory on the edge of a disk containing 
MR-type FQH-liquid should be the same on both sides of 
the plateau. Finally in Section VI, we critically examine 
our results and discuss potential future developments. 
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II. QUASIELECTRON CONSTRUCTIONS 
OVER THE MOORE-READ GROUND STATE 

Here, we briefly review the MR state and its quasihole 
excitations and discuss two existing constructions of trial 
wave functions for quasielectrons and excitons — one us- 
ing the language of CFT, conjectured by Hansson, Her- 
manns, Regnault and Viefers^^'^^, and one conjectured 
by Bernevig and Haldane, who deflne their trial wave 
functions by their vanishing properties'^. Although the 
language used to describe the model states is very dif- 
ferent in the different approaches, the ground state and 
quasihole state wave functions are identical. However, 
each approach has a 'natural' extension towards quasi- 
electrons, leading to distinct, but related model wave 
functions. We will comment on their relation and dif- 
ferences in Sections IV and V. 

In this section, we will focus on the simplest case: the 
bosonic MR state at filling v = 1. In an abuse of lan- 
guage, we still use the words electron and 'quasielectron', 
even though the system is made up of bosons. 



A. Ground state and quasihole excitations 

Let us start by reviewing some important properties 
of the MR ground state and its quasihole excitations. It 
was noted early on in Ref. 38 that Laughlin model states 
as well as their quasihole excitations can be written as 
correlation functions, where the particles are represented 
by CFT operators. Moore and Read generalized this 
approach', and proposed a model wave function based 
on the Ising CFT for the fermionic FQHE state at filling 
2 + 1/2. In the following, we focus on the bosonic version 
of this state — the MR Pfaffian state at filling v = 1: 

^Pf (^1, . . . , z^) = Pf f-^) • Y[{zi - z^). (1) 

V^* ^jJ 

Pf{A) denotes the Pfaffian of the skew-symmetric matrix 
A. It is defined by 

(^) = e^^a(l)a(2)^a(3)a(4)---^a(iV-l)a(iV), (2) 
cr 

where the sum runs over all permutations a of the N 
indices, and e^- is the signature of the permutation. 
The model state (1) is the densest (lowest degree) zero- 

energy state of — the hardcore three-body Hamilto- 
nian (A: = 2), 

ii<---<ik+i i=l 

In particular, the bosonic MR Pfaffian state vanishes as 
the second power of the difference between coordinates 



when three particles come to the same position. More 
precisely, 

N 

^Pf (2:1 = Z2, 2:3, ..,zn) ^ Y[{zi - Zi)^ . (4) 

This vanishing property is a particular case of the more 
general (A:, r) vanishing properties, where the polynomi- 
als vanish as the rth power when A: + 1 particles come to 
the same point 

^^'^^( 51 = ^2 = • • • = ^/c+l, • • • , ^iv) = 
z 

N 

n {z-z,Y^^^':1{z,+,,...,zn). (5) 

i=k-\-l 

Note that for fermionic systems there are ultralocal 
Hamiltonians similar to the T-L^^^^\ which imple- 

ment vanishing properties for the wave functions so that, 
after division of the wave function by a Jastrow factor 
(which is always possible for a fermionic wave function), 
the resulting function still vanishes when A: + 1 of the 
coordinates are equal. 

The quasihole state manifold of the MR state is 
spanned by less dense (higher degree) polynomials that 
satisfy the vanishing conditions (4). These quasihole 
states are in fact ground states (i.e. zero-modes) of the 
three-body Hamiltonian (3), albeit at a higher number of 
ffux quanta A^^ than the Pfaflian state (1). Nevertheless, 
we will call them excitations, because, when perturbing 
away from the model Hamiltonian towards more realis- 
tic Hamiltonians, the degeneracy between these states is 
split and the resulting band should give a good descrip- 
tion of the low-energy sector of the more realistic system. 

One may choose a basis for the space of quasihole states 
which consists of eigenstates of the total angular momen- 
tum operator L and the angular momentum along the 
^-direction, Lz. The quasihole counting, i.e. the number 
of basis states (or equivalently the number of multiplets) 
at each angular momentum, is a fingerprint of the topo- 
logical order of the model state. For the MR Pfaffian 
states, it was explained in Ref. 6 how to calculate the 
number of quasiholes states Af{N, Uqh, lz) for particles 
and Uqh quasiholes with Lz eigenvalue Formulas for 
A/'(7V, Uqh^ lz) for the entire family of Read-Rezayi states 
were obtained using CFT methods in Ref. 25 (see also 
Ref. 40) and by direct counting of polynomials with the 
required vanishing properties and degree restrictions in 
Ref. 26. The number of multiplets at L > can always 
be easily found from the numbers of states at given Lz 
eigenvalue to be M{N^nqh^L) — M{N^nqh^L + 1). To 
give detailed results on Af{N,nqh,lz), let us first define 
the q- binomial [^] by 

"^^1 f , , for a, 6 G N, a > 6 

= <^ {q)a-b{q)b , , — j^g^ 

h otherwise 
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where {q)m = YYJLii^ + Q^)- Following Ref. 25 the gen- 
erating function of J\f{N, Uqh, Iz) can then be written as: 



Nn/A 
l^ = -Nn/4 

N/2 
a=0 



PfafRan ground state can then be written as 



" riqh/2 ■ 
TV -2a 



riqh + a 
a 



(7) 



For the first min[Ar/2, Uqh] angular momentum multiplets 
(counted from the highest) the quasihole counting is iden- 
tical to the edge mode counting that is expected in the 
thermodynamic limit in the disc geometry. This is rea- 
sonable: the sphere and the disc are connected via a 
stereographic projection. The south pole of the sphere is 
mapped to the origin, while the North pole is mapped to 
the edge of the disc. Thus, we expect the state counting 
at high angular momenta on the sphere to correspond to 
the edge counting in the disc geometry. In Section V, we 
will see that the situation is similar for the quasielectron 
state counting; the high angular momentum counting is 
identical to the edge counting, even though finite size 
corrections to the counting appear earlier than for the 
quasihole case. 



B. Quasielectron construction based on CFT 

We now discuss the construction of trial wave functions 
for quasielectrons by Hansson et al.^^ using CFT. Note 
that this construction focused on localized excitations. 
However, by expanding such wave functions in a basis of 
eigenstates of angular momentum, one may always obtain 
a set of trial wave functions for low energy excitations on 
a sphere, which can be compared to the set of trial wave 
functions to be presented in this paper. The Pfaffian 
ground state, Eq. (1), has a natural interpretation as a 
symmetrized wave function of two independent layers of 
Laughlin z/ = 1/2 states: 

^Pf(2;i, ...,zn) = 

( N/2 \ 

S\ n {z,-z,)\z.^,-z.^^f\. (8) 

\i<j=2 j 

Following Cappelli et al.^^'^^, we may now observe that 
this state can be written as a CFT correlator using two 
independent bosonic fields: 0c , which is related to the 
filling fraction and thus to the electric charge, and 0/, 
which distinguishes the two different layers. In the CFT 
we associate an operator V with the electron, given by 



= cos[0Kz)]e^^^^^^ 



(9) 



where V±,{^z^ = exp[zbi0^(z) + i(\)c{z)\ can be interpreted 
as electron operators in layer 1 (+) and layer 2 (— ). The 
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*pf(^i,...,^iv)=(n^(^^-w 



i=i 



= S 



' N/2 N 
j=l j=N/2+l 

(10) 



The homogeneous compensating background charge op- 
erator Obg, see Ref. 3 for details, is needed to render 
the correlation function charge neutral, and thus, non- 
zero. Note that it only contains the field (pc, but not the 
layer field 0^ . This means that non-zero correlation func- 
tions necessarily need to be charge- neutral in <pi. The 
background charge reproduces the correct Gaussian fac- 
tor needed for a valid LLL wave function, but is otherwise 
of no importance for the remainder of this section. 

Using this two-layer description, the non-abelian quasi- 
hole at position rj has to be described by two operators. 



H±{r]) = exp 



(11) 



Without the symmetrization procedure, i7+ and H- 
would be Abelian (Laughlin-type) quasiholes in each 
layer. Because of charge-neutrality, the quasiholes can 
only be inserted in pairs. The non-abelian nature of 
the quasiholes manifests itself in a topological degener- 
acy of 2^~^ for 2n localized quasiholes. In Moore and 
Read's original description using the Ising CFT, this de- 
generacy originates in the two possible fusion channels of 
the CFT operator describing the non- Abelian quasihole. 
In the two-layer description, the degeneracy comes from 
the possible 'distribution' of the 2n quasihole positions 
^1, • • • 5 ^2n over the two layers. This naively over-counts 
the number of quasihole states. However, using tech- 
niques described in Ref. 5, one can show that not all 
of them are linearly independent and that suitable lin- 
ear combinations reproduce the localized quasihole wave 
functions obtained from the Ising description. For in- 
stance, for four quasiholes, there are 2 linearly indepen- 
dent states: 

N 

^Ta' = {H^{m)HMH-{rjs)H_{rj,)l[V{zj)0,,) 

N 

^t^^^ = {H^{m)H_{v2)H^{v3)H-M n n^jW • 

Expanding the localized quasihole states in angular mo- 
mentum eigenstates reproduces the zero-energy multi- 
plets found by diagonalizing the three-body Hamiltonian 
(3). 

In Ref. 31, Hansson et al. introduce quasielectron 
operators which play a similar role in the definition of 
trial wave functions with quasielectrons as the operators 



5 



H± do for wave functions with quasiholes. The guid- 
ing principle used there to construct these quasielectron 
operators is to view the quasielectron as the antiparti- 
cle of the quasihole. However, the operators H^^{r]) 
are obviously not good candidate quasielectron opera- 
tors, because they produce singularities in the electron 
coordinates. Instead, the authors of Ref. 31 constructed 
well-defined, regularized operators V± with the same long 
range properties as H^^{r]); For more details on the reg- 
ularization, see Ref. 32. The operators, V± can be in- 
terpreted as abelian quasielectrons in the =b layers. As 
was the case for quasiholes, explicit symmetrization is 
essential for the non-abelian properties. It was shown 
in Ref. 32 that localized quasielectrons have the same 
topological multiplicity as localized quasiholes, that is, 
for 2n localized quasielectrons, there are 2'^~^ linearly 
independent candidate wave functions. For instance, the 
2-quasielectron candidate wave function is unique and 
given by (see Ref. 32 for the 4- quasielectron candidate): 

N 

/ N/2 \ N/2 

^(mz,+mz^/2+i)/si' Y[{zi - Zj) J| {zi-zjf 



J=2 



2<i<j 



N 



N 



%+i n (^f+i-^j) n i^i-^jf 

J = f +2 / N/2+2<i<j 



(12) 



where we abbreviated dj = dzj. This method can be 
applied for an arbitrary number of quasielectrons and/or 
quasiholes by inserting the appropriate operators in the 
first line of (12). 



Quasielectron construction based on vanishing 
conditions 



Bernevig and Haldane conjectured a quasielectron con- 
struction by imposing vanishing and clustering properties 
on the candidate wave functions^^. Before going into the 
details of their construction, let us review some impor- 
tant background material. The bosonic RR model states, 
as well as their quasihole excitations, are uniquely defined 
by their vanishing properties. In the case of the bosonic 
MR state, the ground state is the lowest degree symmet- 
ric polynomial that vanishes when 3 particles are at the 
same positions. Higher degree polynomials, obeying this 
vanishing property, span the quasihole state manifold. 

The single-particle states in the LLL, 

(l)n{z) = (27^r^!2^)-^^^e-l^l'/^ 

are eigenstates of the angular momentum operator, Lz 
with eigenvalues which are just proportional to powers 
of the complex coordinate z. Hence, there is a basis 



for the many-body states which consists of symmetrized 
monomials mx, where A is a partition of the total angu- 
lar momentum n^. Alternatively, one can label a mono- 
mial by its corresponding occupation number configura- 
tion nx = {nj(A), j = 0, 1, ... , A^^}, where nj{X) is the 
occupation number of the single-particle state with an- 
gular momentum j. A set of partitions may always be 
partially ordered by dominance, denoted by '>'. A parti- 
tion /i dominates another partition u {/n > u) if the latter 
can be obtained by successive squeezing operations on /i. 
Squeezing is a two-particle operation that changes the 
angular momenta of two particles from ji and j2 to j[ 
and such that ji < j[ < ^2 < j2 and total angular 
momentum is conserved. 

It has been realized recently^^'^^, that many frac- 
tional quantum Hall trial wave functions, in particular 
the ground states and quasiholes of the Read Rezayi se- 
ries, can be written as Jack polynomials. This means in 
particular that they have non-zero coefficients only for a 
small subset of the monomials which span the full Hilbert 
space. In fact, for each of these wave functions, there is 
a special partition Aq , called the root partition, such that 
mxQ has non-zero weight and Aq dominates any other 
partition present in the expansion: 



(13) 



Jack polynomials satisfy further properties in addition 
to the fact that they have a nontrivial root partition. In 
particular, the coefficients of the nonzero monomials in 
the expansion of a Jack may be obtained from a recur- 
sion relation. ^^'^^ However, in describing Bernevig and 
Haldane's quasielectron construction, we will use wave 
functions which have a given root partition and addi- 
tional vanishing properties, but which are not necessarily 
Jacks. 

Because the MR ground state is the lowest degree poly- 
nomial that satisfies the vanishing properties (4), and 
inserting quasielectrons necessarily involves lowering the 
total degree, Bernevig and Haldane suggested that the 
quasielectron polynomials are defined by modified van- 
ishing conditions. They conjectured root configurations 
and vanishing conditions for two types of quasielectrons 
which they call Abelian quasielectrons and non- Abelian 
quasielectrons. Abelian quasielectrons carry a full flux 
quantum (that is, in creating Abelian quasielectrons, one 
must lower the electric flux by one quantum per quasi- 
electron), while non- Abelian quasielectrons carry only 
half of a fiux quantum, like the non-Abelian anyonic 
half-fiux quasiholes of the Pfaflaan state. In their paper, 
Ref. 33, Bernevig and Haldane focus on systems with any 
number of Abelian quasielectrons (all localized near the 
same position) and systems with a single non-Abelian 
quasihole and a single non-Abelian quasielectron (these 
are really excitons) . Let us start by reviewing the vanish- 
ing and clustering conditions for the Abelian quasielec- 
trons of the MR state, because they are slightly simpler. 
In the following, the expression 'forming a cluster of n 
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particles' denotes that n particles are at the same posi- 
tions. The requirements on an 5— Abelian quasielectron 
state (which is a state with s Abelian quasielectrons lo- 
calized near the same position) are that it vanishes when 
5 + 1 clusters of 4 particles are formed, and it vanishes 
when one cluster of 2^ + 3 particles is formed as the 
{s + 2)th power of the difference between coordinates. 
In the special case of s = 1 this becomes 

P{zi, Zi, Zi, Zi, Z2, Z2,Z2, Z2, Zq, . . . , Zn) = 

N 

P{zi,...,zi,zq,...,zn) ^ - Zjf (14) 

with the root partition {4 2 2 ... 02}. 

The root configurations and vanishing conditions for 
a non- Abelian quasielectron-quasihole pair for the RR 
Zk states are also given explcitly in Ref. 33. The 
quasielectron-quasihole candidate states for the MR state 
form angular momentum multiplets L = 2, 3, ... , N/2. 
The highest weight states of these multiplets are defined 
by the fact that they have the following root configura- 
tions 

L = N/2; {301111... 111} 
L = 7V/2-l; {30111 1... 102} 



L = 2; {301020. ..202} 



(15) 



and by requiring in addition that they vanish for 2 clus- 
ters of 3 particles and for a single cluster of 4 particles: 



P{zi,Zi,Zi,Z2,Z2,Z2,Z7, . . . , Zn) = 
P{zi, . . . ,^1,2^5, . . . ,ZAr) = 0. 



(16) 



These conditions can be generalized to several quasihole- 
quasielectron pairs and also to an even number of non- 
Abelian quasielectrons^^. The quasielectron-quasihole 
states satisfy a further vanishing property, namely 



P{zi,Zi,Zi,Z2,Z2,Ze, ...,zn) 



{^1-Z2?Y^^,{Z^-Z,)\z2-Z,f 



(17) 



Jain^^. Composite fermions provide an attractive phys- 
ical picture and very successful trial wave functions for 
the most prominent filling fractions in the LLL. For a 
much more in depth review and extensive references, see 
for instance Ref. 46. 

Jain conjectured that states of strongly interacting 
electrons can be understood in terms of states of non- 
interacting or weakly-interacting composite particles, 
called composite fermions. A composite fermion consists 
of a fermion (boson) and an even (odd) number m of 
vortices. When moving in a magnetic field, the attached 
vortices generate a Berry phase that partly cancels the 
Aharonov-Bohm phase. Thus, the CFs behave as if they 
were subject to a reduced magnetic field, B* = B — m^oP 
with p being the density and the magnetic fiux quan- 
tum. The reduced magnetic field B* gives rise to Landau 
like levels, called A levels (ALs) in the following, which 
are separated by an effective CF cyclotron energy hw^. 
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FIG. 1: Sketch of the composite fermion mapping in the case 
of two vortices attached (m = 2 in Eq. (18)). 



Based on this interpretation, Jain proposed a gener- 
alization of the Laughlin wave function, describing the 
FQHE at filling fraction u = effectively as an inte- 

ger quantum Hall (IQH) state of CFs at filling = n: 



This condition, in combination with the second condition 
in (16) also uniquely determines the space of quasihole- 
quasielectron pair wave functions. We compare our trial 
wave functions for excitons and quasielectrons to those 
proposed by Bernevig and Haldane in Section VB. 



III. JAIN COMPOSITE FERMION PICTURE 

A. CF quasielectrons and excitons 

In this section we give a short review of the treat- 
ment of excitations over fractional quantum Hall plateaus 
based on composite fermions (CF), as introduced by 



CF 



_{zi,...,zn) =Vlll 



N 



Y[{Zi - Zj)'^(l)n{zi, Zn) 

(18) 



Here, (/>n(^i, ^at) is the IQH ground state wave func- 
tion with V* = n completely filled ALs. The m Jastrow 



factors ]\f^j{zi — Zj)^ attach m vortices to each particle, 
and Vlll projects the wave function on the lowest Lan- 
dau level (LLL). In the particular case n = 1, the Jain 
state (18) is identical to the Laughlin wave function for 
filling V — 

c> m+l 

An important property of the CF picture is that it 
not only describes the ground states accurately, but also 



gives a very good description of the low-energy excita- 
tions, both neutral (excitons) and charged (quasiholes 
and quasielectrons) . The trial wave functions for excita- 
tions are obtained by creating excitations in the i^* = n 
IQH state of the CFs (again, see e.g. Ref. 46 for a de- 
tailed review). Here, we don't present the method in full 
generality, but rather give some representative examples. 

The first example consists in creating an exciton with 
energy one (in the effective CF cyclotron energy hw'^ unit) 
over the Lauffhlin state with filling fraction v = -^rr 

c) m+1 

(which corresponds to taking n = 1 in (18)). To create 
this excitation, a CF in 0i(^i, ^^at) (see Eq. (18)) is 
removed from the angular momentum / state in the lowest 
A level (LAL) and placed in the angular momentum j 
state in the second A level, thus leaving a hole in the 
LAL (see Fig. 2b). This creates an exciton with total 
eigenvalue j — 

^^^|(zi,..,ZAr) = 



(c) 



i— 







(b) 



(d) 



FIG. 2: (a) The bosonic jy — 1/2 Laughlin state in the com- 
posite fermion picture, (b) An exciton over this state, with 
one CF in the second AL and one hole in the LAL. (c) A 
quasielectron excitation of the Laughlin state: as the number 
of flux quanta is decreased by one unit, one CF has to occupy 
the second AL. (d) A quasihole excitation of the Laughlin 
state: as the number of flux quanta is increased by one unit, 
there is one hole in the LAL. 



with, up to overall normalization. 



(19) 



-i{zi,...,zn,Zi,...,zn) = e*i'--'*" 
xzfzj ...z'rHzi,^,4^^)zl+' 



z^- 



(20) 



where e*^' * '*^ is an antisymmetric tensor, and there is 
an implicit summation over repeated indices. The LLL 
projection can be implemented by putting all z's on the 
left and performing the following replacement in (19)^^: 



dz 



d 



(21) 

'^af~) becomes an oper- 



Therefore, ^i(zi, zat, 2 

ator acting on the product — Zj)^. In a similar 

way, we can generate charged excitations with angular 
momentum A quasielectron state with kinetic energy 
one is obtained by decreasing the number of quantum 
fluxes in 4>i{zi, ...^ zn) by one and placing a CF in the 
angular momentum / orbital in the second A level (see 
Fig. 2c) and finally projecting the wave function to the 
LLL using (21). A quasihole state over the Laughlin state 
is obtained by increasing the number of quantum fluxes 
in 01 (2^1, ...,zn) by one. This creates a hole in the LAL 
with a given angular momentum (see Fig. 2d). 



B. Vanishing properties of composite fermion 
states 

The vanishing properties of the states obtained from 
the CF picture are directly given by the index of the high- 
est A level occupied^^: the bosonic states and the bosonic 



counterpart of the fermionic states (fermionic states di- 
vided by a Vandermonde determinant) constructed from 
n AL satisfy equation (5) with k = n and r = 2. For 
instance, the bosonic Laughlin state and its quasihole 
states (see Figures 2a and 2d) reside completely in the 
LAL; thus, they vanish when two particles are brought 
to the same point. The quasielectron and exciton states 
(see Figures 2b and 2c) involve states in the second A 
level; therefore, they vanish when three particles are 
brought to the same point. Moreover, the CF quasi- 
electron states over a Laughlin state also have vanishing 
properties when multiple clusters are formed simultane- 
ously: the 1-quasielectron states vanish when two clusters 
of two particles are formed^^, whereas the 2-quasielectron 
states vanish when three clusters of two particles are 
formed^^. 



The trial wave functions for quasielectrons based on 
CFT obey the same vanishing conditions for single and 
multiple clusters as above. In fact, we expect the two 
approaches to yield identical vector spaces of trial wave 
functions for quasielectrons. In Ref. 51 this was shown for 
a single localized quasielectron (up to boundary effects 
that are absent in the spherical geometry) . It should hold 
for several quasielectrons as well, when using the exact 
projection to the LLL^^ in the CF construction. For 
neutral excitations of energy > 2, the two constructions 
can in principle differ, as the CF construction involves 
higher ALs. However, for neutral excitations of energy 
two we found that not to be the case. We have verified 
all these assertions numerically for small system sizes. 
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IV. CONSTRUCTING THE ANSATZ WAVE 
FUNCTIONS 

A. Description 

From the composite fermion point of view, we can 
interpret the v = 1/2 Laughhn states in the two-layer 
description of the bosonic MR ground state (8) as Jain 
ground states corresponding to a z/* = 1 IQHE system 
of CFs (consisting of a boson and one attached vortex). 
Therefore we propose to construct the excitations over 
the bosonic MR ground state by creating CF excitations 
over the Laughhn states in each layer of (8). Trial wave 
functions for excitations over the fermionic MR are ob- 
tained by multiplication of the bosonic trial wave func- 
tions with an overall Jastrow factor. 

The physical interpretation of this proposal is compli- 
cated by the explicit symmetrization in (8). A number 
of works addressing this issue have appeared, focusing 
on the fermionic case. It was suggested by T.L. Ho in 
Ref. 52 that the spin-polarized MR states can be obtained 
from the two-layer 331-state^^, whose spatial wave func- 
tion does not include the symmetrization, by introduc- 
ing tunneling between the layers. This proposal turned 
out to be problematic, and it was argued for example in 
Refs. 6,54,55 that tunneling actually drives the system 
to an Abelian phase, although at a special value of the 
parameter that drives the tunneling, a modified version 
of Ho's model (with some three body interaction added) 
does describe a critical theory adjacent to the MR phase. 
Recently, it was suggested in Ref. 55 that a weak MR 
phase can be stabilized around this point, provided that 
one changes the density of the system at the same time 
as introducing the tunneling. While this is clearly an 
important issue, we will not worry too much here about 
the precise physical mechanism that provides the sym- 
metrization (either for bosons or for fermions) , but rather 
simply write down trial wave functions based on the two 
layer picture and subject these to numerical scrutiny. 

Explicitly then, we propose to write a general excita- 
tion over the bosonic MR ground state as: 

S (^^^;''(^i,...,^iv/2) X ^^^;''(^iv/2+i,...,^iv)) (22) 

where ^^^^ and ^^^'^ are excitations (neutral, quasi- 
hole or quasielectron) of the Laughhn state within the CF 
picture. For instance, we can create a two- quasielectron 
state with ^-angular momentum 1^ by taking the ^^p* 
states to be CF quasielectron states with /^^ in layer 1 
and in layer 2 (such that Iz-^ -\- lz2 = h)- 

Notice that in Eq. (22), the projection onto the LLL 
is done separately for and . As mentioned, 

this construction can be extended to the fermionic cases 
by multiplying Eq. (22) by a global Jastrow factor. This 
operation being invertible and L preserving, the number 
of excited states and their angular momentum counting 



obtained in this way for the fermionic states are the same 
as those for bosonic states. 

It is important to note that obtaining the number 
of linearly independent angular momentum multiplets is 
not trivial. As in the case of quasiholes, the trial wave 
functions for quasielectrons that we propose are not all 
linearly independent. The LLL projection in each of the 
layers projects some linear combinations of quasielectron 
states to zero and the symmetrization between the layers 
often introduces further linear dependencies. The same 
considerations apply also to excitons. 

We should also point out that the place where the pro- 
jection onto the LLL is performed is highly relevant. Sev- 
eral others schemes would have been possible. For exam- 
ple, we could have done the projection after the sym- 
metrization. Considering the fermionic states offers even 
more options since projection can be performed before 
or after considering the additional global Jastrow factor. 
Our choice is motivated by the physical picture that the 
bosonic MR state can be seen as two u = 1/2 CF layers, 
both for the ground state and the quasihole excitations. 

An alternative set of trial wave functions for quasi- 
electrons and excitons over the fermionic Pfaffian state, 
presented by Sreejith et al. in Refs. 27,28 can in fact be 
considered as a differing from the trial wave functions we 
propose in the way that the LLL-projection is done. Sree- 
jith et al. define their bilayer composite fermion (BCF) 
wave functions as follows (omitting Gaussian factors), 

a(^^1''^^^{zi,...,zn/2) X ^1''^^^{zn/2+u...,zn)) (23) 

Here the A stands for total antisymmetrization. As in 
our own proposal, the electrons have been split into two 
groups, or layers. However, here each layer is in a z/ = | 
Laughlin (or CF) state, with potentially some quasiholes, 
quasielectrons or excitons, all created according to the 
usual CF construction. There is also repulsion between 
electrons in different layers, but this only induces a sin- 
gle zero in the wave function when two such electrons 
approach. Finally the wave functions is antisymmetrized 
in order to make the electrons all indistinguishable. We 
have taken the number of electrons in each layer equal 
to N/2 in the expression above, but these numbers can 
in principle be different and must be different if the to- 
tal number of electrons is odd (as in Ref. 27). In order 
to compare these wave functions to our own, we divide 
by a Jastrow factor (recall our own wave functions for 
fermions are obtained by multiplying the wave functions 
(22) by a Jastrow factor). This leads to the bosonic BCF 
wave functions 

^BCF,zy=l bosons (^1' ^n) = 

/ ^-JJ^^i^zi,...,ZN/2) '^"JJ'iizN/2+l,--,ZN)\ (24) 

This expression obviously leads to the same groundstate 
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and quasihole states as equation (22). However, for ex- 
citons and quasielectrons, where the LLL-projection in 
the individual layers' wave functions is nontrivial, there 
are differences, because LLL-projection does not com- 
mute with multiplication by a Jastrow factor. Despite 
these differences, one would expect that the two sets of 
wave functions have large overlaps, probably describing 
the same universality class of FQH states. In fact, in nu- 
merical studies of systems with a large number of com- 
posite fermions, one usually does not use the canonical 
LLL-projection, but instead the method introduced by 
Jain and Kamilla in Refs. 56,57. Similarly, for bosonic 
systems with a large number of particles, one typically 
uses the wave functions obtained by first multiplying with 
a Jastrow factor, then projecting using the Jain-Kamilla 
method and then dividing out the Jastrow factor again 
(this practice was introduced in Ref. 58). If we were to 
follow both of these conventions, then the wave functions 
(22) and (24) would become identical. We have in fact 
done our numerical work using the exact LLL projec- 
tion in each layer so that we really study different wave 
functions from those proposed by Sreejith et al. How- 
ever, as a check, we have also done some calculations of 
the number of independent states for small system sizes 
using the alternative LLL-projection which leads to the 
wave functions of Sreejith et al. This gave essentially the 
same results as obtained for our own wave functions, but 
different results from those found by Sreejith et al. More 
detail on this can be found in Section IV C. 



B. Vanishing properties of the trial wave functions 

As the vanishing properties of CF wave functions are 
known (see Section IHB), we can deduce a priori van- 
ishing properties of the two layer states we construct. 
If "^^p^ (^^p^) vanishes when ki {k2) particles are 
brought to the same point, then ^f^f vanishes when a 
cluster of /ci -f /c2 — 1 particles is formed. For example, a 
two quasielectron excitation of the MR state can be built 
from a single quasielectron excitation as depicted in Fig. 
2c) in each layer. Such a state is automatically a zero en- 
ergy state of the Hamiltonian 'H^^ (see Eq. (3)). While 

such states cannot be a zero energy eigenstates of 
since the MR state is the densest zero energy ground state 
of this Hamiltonian, suitable linear combinations of the 
quasielectron states may still be zero energy eigenstates 
of the T-L^^^ Hamiltonian. 

In a similar way, we can deduce the vanishing prop- 
erties when multiple clusters of electrons are formed in 
the states constructed by equation (22) from the cluster 
vanishing properties of the Laughlin quasielectron states. 
When ^^"^'^ and ^'^p' are Laughlin states with 1 quasi- 
electron, the resulting states will vanish when a cluster 
of 4 particles and a cluster of 3 particles are formed and 
when 3 clusters of 3 particles are formed. If ^^p^ and 
are Laughlin states with 2 quasielectron states, the 



resulting states must vanish when a cluster of 4 particles 
and a cluster of 5 particles are formed. This property is 
trivially satisfied since these states already vanish when 
any 5 particles are brought to the same point. How- 
ever, one may also deduce non trivial vanishing proper- 
ties: such states vanish when 2 clusters of 4 and a cluster 
of 3 particles are formed, or when 1 cluster of 4 parti- 
cles and 3 clusters of 3 particles are formed, or when 5 
clusters of 3 particles are formed. 

The idea behind the construction (22) is the same as 
for the CFT construction^^ see Section II B: both are 
inherently two-layer constructions. Due to the similari- 
ties of the two approaches in each layer, we expect that 
they give equivalent descriptions of the low-energy exci- 
tations of the MR ground state. In particular, one can 
show that they obey the same vanishing properties as 
described above, and they yield identical candidate wave 
functions for a single exciton as well as for quasielectrons. 



C. Implementation of state counting and overlap 
calculations 

We use two different methods to generate our trial 
states and compute overlaps. 

The first method calculates the wave functions (22) in 
real space and computes the overlaps using Metropolis in- 
tegration. Real space techniques can be used to compute 
composite fermion wave functions with over 100 parti- 
cles if the LLL projection of Refs. 56,57 is used. How- 
ever, the symmetrization procedure involves a number of 
terms which grows factorially with the number of parti- 
cles N and this limits the reachable size with this method 
to TV = 16. 

In our second method, we first compute Laughlin states 
with the desired excitations in the CF picture us- 
ing the exact method for computing CF wave function 
explained in Ref. 48. Then the states are symmetrized at 
the Fock space level. Using this method, we were able to 
generate bosonic states up to 20 particles. Using Schur 
polynomials that can be generated using recursion formu- 
las from Ref. 59, we can multiply these wave functions 
by a global Jastrow factor at the Fock space level, con- 
verting our bosonic states to fermionic ones. This can 
be done for up to 14 particles (Hilbert spaces of a few 
hundred thousands of independent states). 

A major advantage of the Fock space method over the 
real space method is that all calculations are done at the 
machine precision and results are expressed in the n-body 
basis. In the real space method, overlap calculations are 
done by Monte Carlo integration and suffer from statisti- 
cal errors which are much larger than machine rounding 
errors. A potential advantage of the real space method is 
that it makes multiplication by a Jastrow factor trivial, 
which means fermionic calculations could in principle be 
done for up to 16 particles using this method. 

We now describe our calculation of the number of lin- 
early independent trial wave functions at each value of 
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Lz. Formula (22) provides, for each value of L^, and for 
each number of flux quanta and quasielectrons or exci- 
tons, a set of wave functions z = 1, . . . , g?. These wave 
functions are usually not linearly independent as both 
symmetrization and projection induce linear dependen- 
cies. We want to find a basis of linearly independent 
states x^, i = 1, . . . , (i' < (i for the space spanned by the 
^i. To do this, we compute the overlap matrix, given by 

= i^il^j), (25) 

and we diagonalize it. In the ideal case, the matrix M 
has a number of eigenvalues which are clearly non-zero 
and a number of eigenvalues which equal zero to numeri- 
cal accuracy. This actually occurs in all our calculations 
when we use the Fock space method. Since M is hermi- 
tian, there exists a unitary matrix Q such that QMQ~^ 
is diagonal, that is, ^j^kQij^j^Qki = ^i^ih where 
are the eigenvalues of M. We now define Xi = Qij^ j 
and it follows that (XilXj) = ^ij^i- This implies that 
whenever A^ = 0, we also have x* = 0? giving a linear re- 
lation between the The Xi which belong to nonzero 
eigenvalues, A^ 7^ form the sought after orthogonal ba- 
sis for the vector space spanned by the "^i. Hence, the 
dimension of the space spanned by the trial wave func- 
tions is just the number of nonzero eigenvalues of the 
overlap matrix M. 

The method described above works very well for us 
when we use the Fock space method of evaluating the 
wave functions and overlaps. However, if we use the real 
space method to calculate the overlap matrix, there are 
statistical errors in the matrix which cause its spectrum 
to have a number of spurious small but nonzero eigen- 
values, which should be discarded to obtain the correct 
counting of states. Moreover, it is important to real- 
ize that the dimension of the space spanned by a set of 
trial wave functions is not a quantity that is stable under 
perturbation; random perturbation of a set of linearly 
dependent wave functions will tend to make them all 
linearly independent. Such perturbations may be intro- 
duced for example by making a change in the way that 
lowest Landau projection is implemented. In such cases, 
as long as the spurious eigenvalues of M are small enough, 
one may approximate the original trial wave functions 
very well using only those vectors Xi which correspond to 
the larger eigenvalues. Explicitly, let us first define the 
orthonormal basis vectors Xi = We may expand 

each of the in terms of these, that is, we may write 
= ^jCijXj foi" some coefficients Qj. We then find 
that Cij = (xjl^i) = \^Qji- Since Q is a unitary ma- 
trix, \Qij\ < 1 and hence < Xj. Hence we see that 
it is an excellent approximation, in terms of the quantum 
mechanical inner product, to drop the states correspond- 
ing to the small eigenvalues of M from our description 
of the space spanned by the as long as the sum of 
the M-eigenvalues of the dropped states is much smaller 
than 1. In practice we find in our calculations using the 
real space method that it is always possible to make a cut 



in the spectrum of M which satisfies this condition very 
well, and when this is done we obtain the same counting 
of states as that obtained using the Fock space method. 

Once we have a basis of linearly independent states 
Xi^i = 1, ...c? for the Hilbert space H^^ we can compute 
overlaps between this set of wave functions and another 
set of wave functions obtained, for example, by exact 
diagonalization. To do this, we need to generalize the 
notion of overlap between single states to overlap between 
subspaces. If we have two bases of normalized states 
(pi, i = 1, .., d and 77^, i = 1, d and we are interested to 
know if they generate the same subspace, we can take 
the trace of the operator that projects one of the bases 
into the other, i.e. we can define the (squared) overlap 
between the two subspaces as: 

d 

Overlap = -^^\< (f)i\r]j >\'^ . 
^ i=i 3=1 

(26) 

This quantity (26) is a natural measure for the overlap 
between two different subspaces. In particular, if the 
subspace generated by both bases is the same, the over- 
lap is equal to one, and if they don't generate the same 
subspace it is easy to see that (26) is less than one. 

V. NUMERICAL RESULTS 

In this section, we present the results of our numerical 
calculations. In Section VA we show that our ansatz 
wave functions produce, for the boundary excitations 
on the disc, the same counting formulas for quasiholes 
and quasielectrons. We also conjecture a counting for- 
mula (for any value of Lz) for two and four-quasielectron 
states. In Section VB, we then compute the root parti- 
tions (see Section II C) of our quasielectron states and we 
show that our states share a number of properties with 
the quasielectrons obtained by Bernevig and Haldane in 
Ref. 49. Section VC contains a discussion on the rele- 
vance of unbalanced states, which have different number 
of electrons in the two layers. Finally, in Section VD we 
show how well our construction describes the low energy 
spectrum of both model and realistic Hamiltonians. We 
also show that the subset of our trial states which van- 
ish when a cluster of 4 electrons is formed is particularly 
successful in attaining low variational energies and large 
overlaps with exact wave functions. 



A. Multiplet counting for the trial wave functions 

As discussed in Section IV A both the projection and 
the symmetrization can create linear dependencies be- 
tween the ansatz wave functions (22). We have com- 
puted, for each L^, and for different numbers of particles 
the dimension of the space of linearly independent trial 
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wave functions for the various types of excitations. The 
calculations have been performed up to TV = 20 involving 
Hilbert spaces as large as 6 x 10^. For quasielectrons at 
high values of Lz , these dimensions match perfectly with 
the values predicted by the CFT describing the boundary 
excitations (low energy excitations) of the MR Pfaffian 
phase on the disc (see Section II A) . Therefore our ansatz 



wave functions for quasiparticles show the same topolog- 
ical properties for quasielectrons as for quasiholes. 

In Tables I-III we give the numbers of independent 
states we have found for various numbers quasielectrons 
and excitons, for each value of the total angular momen- 
tum L and for different numbers of electrons. 





N / L 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


n = 2 


14 





1 





1 





1 





1 


























16 


1 





1 





1 





1 





1 
























18 





1 





1 





1 





1 





1 




















n = 4 


14 


2 





4 


1 


4 


2 


3 


1 


2 





1 




















16 


3 





4 


2 


5 


2 


5 


2 


3 


1 


2 





1 
















18 


3 





5 


2 


6 


3 


6 


3 


5 


2 


3 


1 


2 





1 










n = 6 


14 





3 


1 


5 


2 


3 


2 


2 





1 






















16 


3 


1 


6 


4 


8 


4 


7 


3 


4 


2 


2 





1 
















18 





7 


4 


11 


8 


12 


9 


11 


6 


8 


4 


4 


2 


2 





1 










20 


na 


na 


na 


na 


na 


na 


na 


na 


na 


na 


na 


7 


8 


4 


4 


2 


2 





1 


n = 8 


14 


1 





1 





1 
































16 


2 





3 


1 


3 


1 


2 





1 
























18 


4 


1 


6 


4 


8 


4 


7 


3 


4 


2 


2 





1 















TABLE I: Multiplicities of the angular momentum multiplets for n-quasielectron states. The dash symbol indicates that there 
is no state and na stands for non available. The dimensions corresponding to the highest values of L in boldface stabilize when 
we increase the number of particles and matches the dimensions given by the CFT describing the boundary excitations on the 
disc. 



Table I shows the numbers of independent multiplets of 
our trial states with 2 to 8 quasielectrons, for A/" = 14, 
A" = 16 and A" = 18. Note that the dimensions corre- 
sponding to the higher values of L (in boldface) stabilize 
when we increase the number of particles. As explained 
in Section II A, they correspond to the dimensions given 
by the CFT describing the boundary excitations on the 
disc. Using Eq. (7), one may check that the stable mul- 
tiplicities in Table I, are the same as those observed in 
the high L sector of systems with 2 to 8 quasiholes. 



The number of multiplets of single exciton states is just 
the same as the number of single exciton states over one 
of the two layers of composite fermions. This was to be 
expected as the trial wave functions for single excitons 
have an exciton in one composite fermion layer and a 
ground state wave function in the other. We find that a 
system of A" particles has no multiplets at L = and L = 
1 and has a single multiplet at all subsequent L- values up 
to the maximal L- value where a multiplet occurs, which 
is L = A^/2. 

Table II shows the multiplicities corresponding to 
states with a single exciton and two quasielectrons (top) 



N / L 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


12 


2 


3 


6 


7 


8 


6 


7 


5 


3 


2 


1 














14 


1 


5 


7 


10 


10 


11 


9 


9 


7 


5 


3 


2 


1 










16 


2 


5 


10 


11 


14 


14 


14 


12 


12 


9 


7 


5 


3 


2 


1 






12 


2 


3 


6 


7 


9 


9 


10 


8 


7 


5 


3 


2 


1 










14 


1 


5 


7 


10 


11 


13 


12 


13 


11 


9 


7 


5 


3 


2 


1 






16 


3 


5 


9 


12 


15 


15 


17 


16 


16 


14 


12 


9 


7 


5 


3 


2 


1 



TABLE II: Multiplicities of the angular momentum multiplets 
for: a system with one exciton (in the second AL) over a 
two-quasielectron ground state (top) and a system with one 
exciton (in the second AL) over a two-quasihole groundstate 
(bottom). 



and a single exciton and two quasiholes (bottom). As in 
the case of quasielectrons and quasiholes without excitons 
present, we obtain the same numbers for quasielectrons 
and quasiholes at high L, once again confirming the idea 
that the CFT and TQFT describing the quasielectrons 
should be the same as for the quasiholes. We may also 
conjecture a formula for the bold numbers in these tables. 
They equal the integer parts of ^^^"'"^^^^^^ , where p is the 
position of the L-value counting from the highest L for 
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which a multiplet of trial states exists. Alternatively, we 
may characterize these numbers by the generating func- 



N / L 
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2 
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5 


6 


7 


8 


9 


10 


11 


12 


12 


5 


3 


10 


7 


13 


9 


13 


8 


9 


4 


4 


1 


1 


14 


6 


4 


13 


10 


18 


13 


18 


13 


15 


9 


9 


4 


4 


16 


7 


5 


16 


13 


23 


18 


24 


18 


22 


15 


16 


9 


9 



TABLE III: Multiplicities of the angular momentum multi- 
plets for a system containing excitons with effective CF energy 
2. 



the same at all L- values. However, by inspection we note 
that even in these cases, equality of multiplicities be- 
tween quasihole and quasielectron states still holds for 
some L- values beyond the stable ones. 

Some further relations between multiplet countings at 
different numbers of particles and different numbers of 
electrons can be conjectured using particle hole duality 
in the composite fermion A-levels. One may naively con- 
jecture that the same counting should be obtained for a 
system with particles and n quasielectrons as for a sys- 
tem with 2N — 3n + 4 particles and 2{N/2 — n + 2) quasi- 
electrons (this is true before LLL-projection and sym- 
metrization) . We observe in our data that this holds for 
n = 2 and n = 4, but not beyond 4 quasielectrons. 



Table III shows the multiplicities corresponding to a sys- 
tem with up to two excitons. Observe in the table that 
the stable multiplicities occur in pairs and that they are 
simply the squares of the integers (1, 4, 9, 16 etc.). Dou- 
ble excitons can in principle occur in three different ways. 
One may excite a particle to the first unoccupied AL in 
each of the composite fermion layers, or one may excite 
two particles in a single layer, or one may excite a sin- 
gle particle up to the second unoccupied AL. We have 
considered all of these possibilities, as they all occur at 
the same naive composite fermion energy. However, it 
is worth noting that, to obtain the entire set of linearly 
independent trial wave functions, one does not actually 
need to include the states with two excitons in the same 
layer or with particles excited to the second unoccupied 
AL. We checked numerically that the states constructed 
in this way are already contained in the space of trial 
wave functions with a single exciton in each layer. 

We have already conjectured counting formulae for the 
stable numbers of multiplets of quasielectron and exci- 
ton states at high L, relating quasihole and quasielec- 
tron states. However, for 2 and 4- quasielectron states, 
we may go further and conjecture counting formulae for 
any L- value. Observe that, before projection and sym- 
metrization, the number of independent n-quasielectron 
states (considering second AL quasielectrons only) in a 
system with N particles, is the same as the number of 
independent n-quasihole states in a system with par- 
ticles, where N' is given by: 

7V' = Ar-2*n-h4. (27) 

Of course the naive counting of states is modified as 
the projection and symmetrization operations introduce 
linear dependencies between the ansatz wave functions. 
Nevertheless, for all cases we have checked, it turns out 
that after projection and symmetrization, these count- 
ings, though modified by projection and symmetrization, 
are still equal for quasiholes and quasielectrons, for all L, 
in the particular cases of 2 and 4-quasihole/quasielectron 
states (we checked this up to N = 18). We conjecture 
that this equality holds for all N. For higher numbers 
of quasiholes/quasielectrons, the countings are no longer 



B. Vanishing properties and comparison to 
Bernevig-Haldane construction 

Now we come back to the Bernevig-Haldane quasi- 
electron and exciton states introduced in Section II C 
and we will compare them with our ansatz wave func- 
tions. The easiest comparison to make is between our 1- 
exciton states and those of Bernevig and Haldane which 
have a single non-Abelian quasielectron and a single non- 
Abelian quasihole. We have checked that our single exci- 
ton states have the same root partitions as the Bernevig- 
Haldane excitons (given explicitly in formula (15)). We 
also find that they satisfy the same vanishing proper- 
ties (16) and (17). Hence, for single excitons, our trial 
wave functions are in fact the same as those proposed by 
Bernevig and Haldane. 

For wave functions with multiple excitons or with only 
quasielectrons, it is a bit more complicated to make a 
comparison between our trial wave functions and the ones 
proposed by Bernevig and Haldane, if only because the 
root partitions for such states are not given explicitly in 
Ref. 33. However, it is clear that, in these more general 
cases there can be some mismatch between the two con- 
stuctions. For example, let us compare our construction 
to Bernevig and Haldane's construction for non-Abelian 
quasielectrons. 

As discussed in Section IV B, our two-quasielectron 
states vanish when 5 particles cluster together. They 
also vanish when 3 clusters of 3 particles are formed and 
when 2 clusters, one of 4 particles and one of 3 particles 
are formed. These vanishing properties are also satis- 
fied by our 2-exciton states. However, the non-Abelian 
quasielectron states of Ref. 33 vanish already when a 
single cluster of 4 particles is formed and also when 2 
clusters of 3 particles are formed. It is easy to check di- 
rectly that not all of our trial wave functions satisfy these 
stronger vanishing properties, so we obtain a mismatch 
with Bernevig and Haldane's construction. 

It is interesting to look at subspaces of our space of 
trial wave functions which do satisfy stronger vanishing 
properties, such as those required by Bernevig and Hal- 
dane. One motivation for this is that it may be (naively) 



expected that trial wave functions which vanish aheady 
when 4 particles cluster together may have lower varia- 
tional energy (for repulsive potentials) than wave func- 
tions which don't vanish until a cluster of 5 particles is 
formed. This indeed turns out to be the case (see Sec- 
tion VD). We have therefore studied in some detail the 
subspaces of our spaces of trial wave functions which con- 
sists of those states which vanish when a cluster of 4 par- 
ticles is formed. In other words, these are the space of 
zero modes of the T-L^^^ Hamiltonian inside our spaces of 
trial wave functions. We give a sample of the results for 
the counting of such states in Table IV. By comparison 
to Table I we see immediately that, indeed, not all our 
trial states are zero modes of 'H^^. In particular, the 
trial states at the highest angular momenta seem never 
to vanish when 4 particles cluster together. 





N / L 
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14 
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1 






n = 4: 


14 
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2 
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2 
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TABLE IV: Multiplicities of angular momentum multiplets 
for n-quasielectron states which are also zero modes of 'H^^ . 



We find similar results to those presented in Table IV for 
the counting of multiplets of 2-exciton states which are 
also zero modes of H^'' . Again, not all our trial states for 

double excitons are zero modes of 'H^^ and in particular, 
the trial states at the two highest angular momenta never 
vanish when 4 particles cluster together. For the case 
of double excitons, we can give an exact description of 
the numerical results for multiplet counting by a formula 
which relates the number of multiplets which are zero 
modes of to the total number of multiplets. We 
observe that the total number of multiplets equals the 
number of multiplets of zero modes for L = 0, 1, 2, 3 
and then from L = 4 upward, the number of multiplets 
of 'H^^ zero modes is lower than the total number of 
multiplets, by min{[L/2- 1], [(iV-L)/2-|-l]}, where min 
denotes the minimum and the square brackets denote the 
integer part. 

We also calculated the root partitions for the two- 
quasielectron states which are also zero modes of 'H^^ 
For N = 16, these are given in Table V. These 
root configurations correspond to the ones predicted by 
Bernevig and Haldane's construction^^ for states with 2 
non Abelian quasielectrons. 

C. Unbalanced states 

So far, we have considered only balanced states, that 
is, states with the same number of electrons in both CF 
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highest root configuration 





30111111111103 
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21111111111103 
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20211111111103 
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20202111111103 
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20202021111103 


5 


20202020211103 


6 


20202020202103 



TABLE V: Highest root configurations of the zero energy 
states of 'H^^ in each Lz sector of our space of 2-quasielectron 
trial states at A/" = 16 

layers. However, there is no strong reason for this re- 
striction a priori. In fact for odd numbers of electrons, 
unbalanced states with an odd number of electrons in 
one layer and an even number in the other layer are un- 
avoidable, and as shown in Ref. 27 these indeed give a 
good description of the lowest band of states for an odd 
number of electrons in the MR Pfaffian phase. In this 
paper, we are looking at systems where the total num- 
ber of electrons is even, and in this case, we find that, 
at least in the low energy sector of the theory (defined 
with respect to the CF energy), unbalanced states don't 
introduce new physics and can be ignored. The rest of 
this section is devoted to explaining this in some detail 
and can be skipped on first reading. 

We will deal with unbalanced quasihole states in the 
bosonic case, but the same results hold for fermions. 
Consider a 2n-quasihole state with N = Ni + N2 par- 
ticles where Ni and A^2 are the number of particles in 
the first and second layers respectively and suppose that 
Ni > N2. We will require that both layers feel the same 
magnetic field. This requirement can be physically mo- 
tivated by the idea that the Pfaffian phase can really 
be viewed as some kind of modified two- layer system. 
The largest exponent of the coordinates in the first layer, 
A/'^i, must then be equal to the corresponding exponent 
N(jj2 in the second layer. Note also that these two expo- 
nents must equal the total number of fiuxes in the system, 
N(f, = N — 2 -\- n. Hence, we have 

N^i = N-2^n = 2Ni-2 + {N2-Ni+n) 
N^2 = N-2^n = 2N2-2^{Ni-N2^ n). (28) 

From this equation, it is clear that the second layer cor- 
responds to a = 1/2 Laughlin state with Ni — N2 -\- n 
quasiholes. For the first layer we have two different cases 
(depending on the values of Ni and A^2) corresponding to 
a z/ = 1/2 Laughlin state with either Ni — N2-\-n quasi- 
holes (if A^2 — ^1 + ^ > ) or A^i — A^2 + ^ quasielectrons 
(if N2 — Ni-\-n<0). In the first case the two layers of the 
trial wave function both contain z^ = 1/2 Laughlin quasi- 
hole wave functions. These 2n-quasihole states have the 
property that they vanish when three particles are at the 
same position (see Section II C). Therefore because the 
balanced 2n-quasiholes already span the complete set of 
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wave functions with this vanishing property we can con- 
clude that the unbalanced quasiholes are linear combina- 
tions of the balanced ones and we don't need to consider 
them in our construction. 

The case where one layer has quasiholes and the other 
quasielectrons is more interesting. In particular, in cases 
with a single quasielectron in one of the layers, these 
unbalanced states have the property that they vanish 
whenever a cluster of 4 particles is formed. This vanish- 
ing property is different from the usual vanishing prop- 
erty for balanced quasihole states, which already van- 
ish when 3 particles form a cluster. We have verified 
that these unbalanced states are linearly independent of 
the balanced wave functions and therefore they are new 
quasihole states. However, we have two reasons to be- 
lieve that these unbalanced quasiholes are not relevant 
to the lowest energy sector of quasihole states. First of 
all, we can use the composite fermions' kinetic energy, 
read off from the number of particles in each AL, as a 
guide. The new states include higher AL contributions, 
while the balanced quasiholes have composite fermions 
only in the lowest AL, and so, at least naively, we can 
focus on balanced quasiholes in any low energy descrip- 
tion. Secondly, the fact that the new quasihole wave 
functions vanish only when 4 particle positions coincide 
will probably mean they have higher variational ener- 
gies for realistic repulsive interactions when compared to 
balanced quasihole states which vanish already when 3 
particle positions coincide. 

For unbalanced quasielectron states it can be shown, 
in a similar way as we did for unbalanced quasiholes, 
that depending on the values of A^i and we can have 
either quasiholes or quasielectrons in the second layer, 
while in the first layer we always have quasielectrons. It 
is easy to see that, in contrast to the situation for un- 
balanced quasiholes, the unbalanced quasielectrons can, 
for suitable values of A^i, A^2 and n, produce states with 
the same AL energy as the balanced ones. However, it 
is also easy to see, that this only starts to occur when 
the overall number of quasielectrons is at least 4. In the 
4-quasielectron sector we have compared balanced and 
unbalanced 4-quasielectron states with the same AL en- 
ergy and we have found that the unbalanced states are 
linear combinations of the balanced ones. Therefore they 
can be omitted from our construction. We have similarly 
investigated a number of other cases where unbalanced 
quasielectron states with the same AL energy as the bal- 
anced quasielectron states exist, and in all cases consid- 
ered, the unbalanced states were already contained in the 
space spanned by the balanced states. 

Finally for the exciton sector (where n = 0) it is easy 
to see that unbalanced excitons correspond to states with 
Ni — N2 quasielectrons in the first layer and with the 
same number of quasiholes in the second layer, in ad- 
dition to any excitons that may exist within the layers 
(this was also noted in Ref. 27). It can be shown that, 
also in this case, for certain values of A^i and N2 it is 
possible to construct low energy unbalanced states with 



the same A level energy as the balanced excitons studied 
before. For single excitons this does not happen, as long 
as N is even, but for double excitons, with AL energy 2, 
there are unbalanced states at the same AL energy with 
A^i — 7V2 = 2: 2 quasielectrons in the first layer and 2 
quasiholes in the second layer. However, we checked that 
these unbalanced excitons are linear combinations of the 
balanced two- exciton states and once again we do not 
need to include them in our construction. 



D. Comparison of Spectra and Overlaps 

To test our construction, we now compare the exact 
spectra of realistic and model Hamiltonians in the LLL 
with the spectra of these same Hamiltonians in the spaces 
spanned by our trial wave functions. We do this both for 
bosons and for fermions. For the bosonic states, we use 
the H^^^ Hamiltonian, for which the MR state and its 
quasihole states are zero-energy eigenstates, as our model 
Hamiltonian. We use the Hamiltonian as our real- 
istic Hamiltonian. This is well justified, since in most 
experiments with ultra-cold bosonic gases, the relevant 
interaction is s-wave scattering^^, which is modeled well 
by this potential. For the fermionic states, the realistic 
interaction is the Coulomb interaction in the second Lan- 
dau level. We considered this interaction with the first 
relevant pseudopotential Vi (which describes the shortest 
range part of the interaction^^) shifted by SVi = 0.035, 
so that the overlap between the MR ground state and 
the exact ground state is maximal. Since the fermionic 
MR state and its quasiholes states are zero-energy eigen- 
vectors of the 3-body hollow core interaction^^ (the 

analogue of for fermionic systems), we take this as 
our model Hamiltonian for the fermionic states. Studies 
of spectra and overlaps comparable to the one presented 
here can be found in Ref. 23 for the bosonic quasiholes 
states and in Ref. 22 for the fermionic ones. A recent pa- 
per by Sreejith et al, Ref. 28, also presents a number of 
results on spectra and overlaps closely related to those 
presented here, for a slightly different family of wave 
functions (see Section IV A for details). 

In Fig. 3, we show spectra for systems with 2 quasi- 
electrons. The top panels of Fig. 3 show the spectra 
of 'Hg^ and n^^\ in the fuh Hilbert space and in the 
space spanned by our trial states for 2 quasielectrons, for 
N = 16 bosons at A^^ = 13. The lower panels of Fig. 3 

show the spectra of H^^ and of the second LL Coulomb 
Hamiltonian with slightly modified Vi pseudopotential, 
again in the full Hilbert space and in the space spanned 
by our trial states for 2 quasielectrons, now for = 14 
fermions at N^p = 23. All panels also show the spectrum 
of the relevant Hamiltonians in the space of trial states 
for 2 quasielectrons which have the additional property 
that they are zero modes of H^^^ (for bosons) or of 
(for fermions). 

In all cases, the low- lying part of the spectrum in the 
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0.927 
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0.891 


0.546 
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0.480 
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0.497 
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0.984 


0.798 
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0.935 


0.672 


0.520 


0.478 


0.964 


0.809 


0.985 


0.760 
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0.378 


0.348 
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0.984 


0.798 
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0.312 


0.410 


0.412 
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TABLE VI: Overlap between the space spanned by quasielectron states for N = 16 and the corresponding low energy space 
with respect to the l-L^ and Ti^ Hamiltonian in each L sector. '4body' indicates that only states are taken into account that 
vanish when 4 particles are at the same position. For 2 quasielectron states, the only state that does not have this vanishing 
properties is the state at L = 8. We can notice that overlaps with this state are much smaller than the other ones. A dash 
means that there is no state for the corresponding L value. 



space of trial wave functions is similar to the low ly- 
ing part of the full spectrum. However, for bosons, the 
quasielectron states at the highest angular momentum 
(in this case L = 8) have anomalously large energy and 
do not obviously match anything in the exact spectra. 
For fermions the highest angular momentum states (now 
at L = 7) also have relatively high energy, though the dif- 
ference with the other trial states is not as pronounced 
as in the case of bosons. These high angular momen- 
tum states are also the only trial wave functions which 
do not vanish when 4 particles are brought to the same 
point. More precisely, in the case of fermionic trial wave 
functions, the corresponding bosonic trial wave function 
does not vanish. Thus, it seems that the 4-body vanish- 
ing property satisfied by the trial wave functions at lower 
L may play a role in obtaining good agreement between 
trial states and low-energy states, as predicted by Hal- 
dane and Bernevig.^^. Note that all trial wave functions 
vanish when 5 particle positions coincide. We may also 
note that for bosons, the low-energy part of the spec- 
trum of the model Hamiltonian 'H^^ is obviously better 

reproduced than that of the realistic Hamiltonian. 
For fermions, the quality of approximation of the realistic 
and model Hamiltonians is comparable. 

We have also calculated overlaps between the trial 
states and the corresponding low lying states of the 



Hamiltonian. The results for bosonic states can be found 
in Table VI and the results for fermions in Table VIL 
In both cases, the results for 2-quasielectron states are 
in the first two columns. When we compare our bosonic 
trial wave functions for 2 quasielectrons to the lowest 
eigenstates, the overlaps are all over 0.93, except 
for the 1/ = 8 state for which the overlap is 0.3. The 
overlaps with the lowest 'H^^ eigenstates are lower across 
the board, except at the highest angular momentum, but 
there the overlap is still low at 0.4. This result is consis- 
tent with the fact that the spectrum of 'H^^ is reproduced 

quite a bit worse than that of (see figure 3). For 
fermions, the overlaps of the trial wave functions with 
the spectrum of the model Hamiltonian T-L^^ are not as 
high as the overlaps for bosons, with the lowest overlap 
equal to 0.83 if we exclude the L = 7 state which is not 
a zero mode of However, in the fermionic case, 

the overlaps for the realistic Coulomb Hamiltonian are 
comparable to those for the model Hamiltonian. Also, 
the potentially anomalous state at the highest L-value 
(here L = 7) has considerably better overlap in the case 
of fermions, reaching 0.62 for the Coulomb Hamiltonian. 
Nevertheless it still has much lower overlap than the trial 
states at lower L- values. 



We similarly investigated systems with 4 quasielectrons. 
In Fig. 4, we show the spectra for these systems. The top 
panels again show the spectra of 7/^^ and , in the full 
Hilbert space and in the space spanned by our trial states 



for 4 quasielectrons, for N = 16 bosons at N^p = 12. The 

lower panels show the spectra of Hp^ and of the sec- 
ond LL Coulomb Hamiltonian for TV = 14 fermions at 
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0.955 


0.973 


- 


- 


1 


0.903 


0.792 


- 


- 


- 


- 


- 


- 


2 


- 


- 


0.791 


0.757 
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0.568 


0.911 
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0.569 
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TABLE VII: Overlap between the space spanned by fermionic quasielectron states for = 14 and the corresponding space of 
low energy states with respect to Tip'' and Tic, the Coulomb interaction in the second LL Hamiltonian -\-5Vi — 0.035 in each 
L sector. A dash means that there is no state for the corresponding L value. 



N(j) = 23. In both cases, the number of trial states gen- 
erated is much higher than for 2 quasielectrons. Only a 
few of these states vanish when a cluster of 4 particles is 
formed. As in the case of 2 quasielectrons, these states 
belong to the low energy part of the spectra. Overlaps be- 
tween fixed L subspaces spanned by 4- quasielectron trial 
states and subspaces of lowest energy states of the Hamil- 
tonians at the same L-value are given in the third and 
fourth columns of Tables VI and VII. Overlaps between 
the subspaces of 4-body vanishing states and subspaces 
of lowest energy states of the Hamiltonians are shown in 
the fifth and sixth columns of these tables. The overlaps 
for the 4-body vanishing states are considerably higher 
than those for the full set of trial states and indeed it is 
also easy to see from the spectra that the 4-body vanish- 
ing states are among the lowest energy trial states and 
certainly give a very economical description of the lowest- 
energy part of the full spectrum. For bosons there even 
appears to be a low energy 'band' of states at L = 0, 
L = 2, L = 4 and L = 6 which can be described by a 
subset of the trial wave functions with the 4-body van- 
ishing property. On the other hand it is obvious from the 
spectra, especially from those for fermions, that the full 
set of trial wave functions does give a reasonable descrip- 
tion of a much larger part of the low energy spectrum 
than the 4-body vanishing states. 

We have examined our construction of exciton states 
by the same means. We consider the two lowest effective 
cyclotron energies here, i.e. exciton states with effective 
composite fermion energy up to 2 in units of the effective 
cyclotron energy of the composite fermions. The spectra 
of the different Hamiltonians in the spaces spanned by 
the 1-exciton and 2-exciton trial states and in the full 
Hilbert space are shown in Fig. 5 for bosons and in Fig. 6 
for fermions. Overlaps between single exciton states and 
low energy states in the full Hilbert space are given in 
Table VI for bosons and in Table VII for fermions. 

Single excitons (with effective energy 1) are obtained 



when one of the liquids is in the Laughlin state while the 
second is in a state excited over the Laughlin state with 
one exciton. These states naturally vanish when 4 parti- 
cles are brought to the same point (see Section IV B). In 

the case of the model Hamiltonians l-L^^ and l-L^ ^ the 
magneto-roton-like branch is remarkably well reproduced 
by these 1-exciton states, both in the spectra and in the 
overlaps. For the realistic Hamiltonians the performance 
of the 1-exciton trial wave functions is less impressive, 
but still quite reasonable. 

As discussed in Section V A, all 2-exciton states can be 
generated by considering only the case when both liquids 
are in a 1-exciton state. In this case, the a priori van- 
ishing property is a 5-body cancellation. However, most 
of the states also vanish when 4 particle are brought to 
the same point (see Section VA for details). The dif- 
ferent spectra in the space of 2-exciton states with the 
4-body vanishing property are also shown in Figs. 5 and 
6. Due to the large number of trial states generated and 
the large number of exact low energy states of the various 
Hamiltonians involved, it is difficult to make a meaning- 
ful quantitative comparison of the different spectra, and 
so we have not listed overlaps for 2-exciton states. How- 
ever, we can see from the spectra that our construction 
behaves well with respect to the energy of model and 
realistic Hamiltonians, in the sense that, as we consider 
more trial states we manage to describe more of the lower 
energy part of the spectra. 



VI. DISCUSSION 

We have proposed trial wave functions for charged and 
neutral excitations of the MR Pfafiian phase, tested these 
numerically and compared our construction to existing 
proposals. 

Despite what seems like a very different method of con- 
struction, we find that our wave functions are precisely 
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the angular momentum eigenstates which appear in the 
trial wave functions for localized quasiparticles and exci- 
tons proposed by Hansson et al.^^. On comparing with 
the trial states based on vanishing properties proposed by 
Haldane and Bernevig^^, we find that our single exciton 
wave functions coincide with theirs, but our quasielec- 
tron wave functions do not. Finally, our construction is 
very similar to a CF based construction by Sreejith et 
al.^^ and we show that the difference can be viewed as a 
change in the LLL projection. 

As one of our main results, we find that the counting 
of large L multiplets of independent trial wave functions 
in our construction is the same for quasholes and quasi- 
elect rons, supporting the idea that quasiholes and quasi- 
electrons are described by the same CFT and TQFT. 
This result is in disagreement with the claim of Sreejith 
et al. in Ref. 28 that the counting of quasielectron states 
is different from that of quasiholes. While the trial wave 
functions of Ref. 28 differ from ours in the details of the 
LLL-projection, this is not the explanation of the dis- 
agreement, as we have checked for small systems that, 
with the method detailed in Section IV C, we find the 
same quasielectron multiplet countings for both types of 
trial states. 

We also tested how well our trial wave functions repro- 
duce the spectra of idealized 3-body Hamiltonians and 
of more realistic Hamiltonians, both for bosons and for 
fermions. We find that the low lying parts of the en- 
ergy spectra for the 3-body Hamiltonians are reproduced 
very well when diagonalizing these Hamiltonians in the 
spaces of trial states for quasielectrons and excitons. This 
is also reflected in high overlaps between the low lying 
eigenstates in the spectra, especially for 2-quasielectron 
states and single exciton states. For bosons, the trial 
wave functions reproduce the spectra quite a bit less well 
in the case of the more realistic Tv^ Hamiltonian, though 
the agreement between the spectra in the full space and 
trial spaces and the overlaps are still reasonable. For 
fermions, it appears that the agreement between our trial 
wave functions and the eigenstates of the second Landau 
level Coulomb Hamiltonian (with slightly shifted 5Vi to 
stabilize a MR Pfaffian phase) is roughly equally good 
as the agreement between our trial states and the eigen- 
states of the model 3-body Hamiltonian. 

It would be of interest to perform a detailed compari- 
son between our results for the spectra and overlaps of the 
fermionic wave functions and those in Ref. 28, but this 
is not completely straightforward. Sreejith et al. use the 
Coulomb Hamiltonian without the shift in the pseudopo- 
tential Vi which we use to stabilize the Pfafiian phase. 
Also, the overlaps presented in Ref. 28 are overlaps be- 
tween the lowest energy state of the Hamiltonian in the 
space of trial wave functions at a given value of L and the 
lowest energy state of the Hamiltonian in the full Hilbert 
space at the same value of L. We present instead the 
overlap (26) between the full space of trial wave functions 
at a given L-value and a space of low energy eigenstates 
of the Hamiltonian which has the same dimension. We 



think that our approach is more consistent to probe the 
accuracy of an excitation manifold construction. Never- 
theless, we checked, up to the definition of the overlap 
(we use I (^^1^^) p instead of | (^^1^^) |), that our re- 
sults are in agreement with those of Ref. 28 for the 2 and 

(3) 

4 quasielectrons case with the Hamiltonian. 

While our trial wave functions do reproduce the low ly- 
ing parts of the spectra with at least reasonable success 
for each of the Hamiltonians considered, it is also clear 
that many of the trial states have high variational ener- 
gies and could be dispensed with. In fact, we observe that 
there is a subspace of our space of trial wave functions 
which typically have much better variational energies and 
overlaps compared to the other trial wave functions, espe- 
cially in the case of bosons. This is the space of trial wave 
functions with 4-body vanishing properties, or more pre- 
cisely, the space of zero modes of T-L^^ for bosons and the 

space of zero modes of 'H^^ for fermions. This observa- 
tion suggests that it is a good idea to construct trial wave 
functions based on their vanishing properties, a technique 
already employed by Bernevig and Haldane in Ref. 33. It 
is important to note however, that one must employ more 
complicated vanishing requirements than those given by 
Bernevig and Haldane if one wishes to describe arbitrar- 
ily large numbers of quasielectrons. One reason for this 
is that, as one decreases the number of flux quanta TV^ 
in the system (at fixed A^), from the MR ground state's 
flux of Nff) = N — 2 (for bosons), one eventually reaches 
values of AT^ where no states exist which vanish when 
4 particle positions coincide. This happens for bosons 
when A^0 < |A/' — 2, the flux of the bosonic = 3 Read- 
Rezayi state. One might argue that, in a real system. A" 
is very large and one would never have to worry about 
systems with this many quasielectrons. 

However, in the small systems that are used in numer- 
ical studies, the number of 4-body vanishing states can 
start decreasing already at a relatively modest number 
of quasielectrons (especially if further vanishing require- 
ments are imposed) and thus one cannot hope to extract 
the correct large A" limit of the counting of multiplets of 
quasielectrons from numerical calculations on such states. 
Also, it must be noted that for fermions the difference 
in overlaps and energies between states with the 4-body 
vanishing property and other low-lying trial states is less 
pronounced than for bosons. It seems that the 4-body 
vanishing property does guarantee good variational en- 
ergy, but absence of this property does not mean that 
the energy or overlap will necessarily be bad. 

We will now make some remarks on future directions 
for this line of research. All the work done here can be 
repeated for systems with an odd number of electrons 
(necessarily with different numbers of electrons in each 
layer). No results on multiplet counting have so far been 
published for such systems. For fermions, overlaps and 
spectra for the lowest lying neutral states will likely be 
similar to those presented in Ref. 27. 

All constructions presented here can be generalized 
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straigthforwardly to the Bonderson-Slingerland hierar- 
chy states^^'^^, which potentiahy describe ah observed 
fining fractions for fermions in the second Landau level. 
The ground state wave functions of the Bonderson- 
Slingerland states contain the bosonic Pfaffian wave func- 
tion as a factor; hence quasiparticles and excitons on this 
Pfaffian factor provide a natural set of trial wave func- 
tions for excitations of the Bonderson-Slingerland states. 
Other states, to which our construction generalizes nat- 
urally, are the non-Abelian condensate state of Ref.^^. 
These can, in fact, be interpreted as double layer states 
with a Jain CF state in each layer, even though their 
construction is via CFT. The BCF states proposed in 
Ref. 27 are related to the states in Ref.^^ by a change of 
LLL-projection, in the same way that their pfaffian trial 
wave functions relate to the trial states presented here. 

We therefore expect, in line with the predictions of 
Ref.^^, that they will have quasiparticle multiplet count- 
ings at large L described by parafermionic CFTs. 

Similarly, by working with k composite fermion lay- 
ers instead of just 2 (still with symmetrisation over all 
layers), the construction presented here can be used to 
produce candidate quasiparticle and exciton wave func- 
tions for the entire series of Read-Rezayi states. When 
A: > 2, there is actually still an open question with regard 
to quasiholes in this case; while it is very plausible^^, it 
has not been strictly proven that the quasihole wave func- 
tions produced using the /c-layer picture of these states 
and those that come from the description of the states 
in terms of vanishing properties are equivalent. That is, 
the layered wave functions satisfy the required vanishing 
properties, but is has not been proven that the angular 
momentum eigenstates which occur in their expansion 
span the full space of states with these vanishing proper- 
ties. 

One may also define a version of the bilayer composite 
fermion states of Ref. 28 which utilizes a LLL-projection 
more similar to the one used in the current work. The 
bosonic versions of these wave functions would simply be 
the symmetrized product of two composite fermion wave 
functions, where each of the composite fermion layers is 
at an effective filling i/* = n, with n > 1 integer (the 
case n = 1 with single flux attachment was studied in 
this paper) . It would be of particular interest to find the 
counting of quasiparticle multiplets for these states. 

Another important direction for future work is to test 
our trial states against the spectra of more physically 
realistic Hamiltonians, often necessarily in larger Hilbert 
spaces. Here, one may think for instance of including 
the effects of the electron spin, Landau level mixing and 
subbands in finite thickness quantum wells. 
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